Angular Schmidt Modes in Spontaneous Parametric Down-Conversion 
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We report a proof-of-principle experiment demonstrating that appropriately chosen set of Hermite- 
Gaussian modes constitutes a Schmidt decomposition for transverse momentum states of biphotons 
generated in the process of spontaneous parametric down conversion. We experimentally realize 
projective measurements in Schmidt basis and observe correlations between appropriate pairs of 
modes. We perform tomographical state reconstruction in the Schmidt basis, by direct measurement 
of single-photon density matrix eigenvalues. 
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Introduction. Characterizing entanglement in 
high-dimensional quantum systems is a hot topic 
in quantum information science. One of the less 
studied subjects is entanglement in the extreme 
case of infinite-dimensional by-partite states. It is, 
nevertheless, of much interest, since such states are 
usual for two-particle scattering processes, such as 
spontaneous parametric down conversion (SPDC), 
being one of the most experimentally attractive 
source of entangled photons. An approach to 
quantifying entanglement in this case was pro- 
posed by Law and Eberly It is based on 
Schmidt decomposition - a representation of a bi- 
partite state vector as a sum of factorized terms: 
1*12) = Z^feV^fc |wfc> \vk), with \u k ),\v k ) being 
eigenvectors of reduced single-particle density ma- 
trices (so called Schmidt modes) and X k - corre- 
sponding eigenvalues. This decomposition has sev- 
eral remarkable features: 

• Schmidt modes, being pure single-particle 
states, form a complete and orthogonal basis; 

• decomposition has a discrete single-sum 
form, i.e. every Schmidt mode is correlated 
with exactly one counterpart; 

• number of significant eigenvalues, defined 
as if = 1/ ^2 k Xl , quantifies entanglement 
present in the system. 

In this Letter we experimentally address the phys- 
ical properties corresponding to the listed features 
of Schmidt decomposition. Specifically, we choose 
spatial degrees of freedom of biphoton field as an 
experimentally convenient example of an infinite- 
dimensional system. 

Angular spectrum of biphotons is well recog- 
nized as an attractive object to study multidimen- 
sional entanglement and its properties were ex- 
perimentally investigated in numerous works dur- 
ing the last decade. These were focusing on two 
main alternatives: EPR-like correlations in trans- 
verse momentum [2|45| and entanglement in orbital 
angular momentum [6t4l7j| . Both approaches may 
be addressed on the same grounds: transverse en- 
tanglement manifests itself in non-classical correla- 
tions between coherent spatial modes - plain waves 



in the first case, and arbitrarily chosen Laguerre- 
Gaussian modes in the second one. None of these 
choices is perfect in the following sense: the cor- 
relations between modes are non-ideal, meaning 
that a single mode is correlated to multiple coun- 
terparts. For EPR case the number of correlated 
modes is limited from below by finite angular di- 
vergence of the pump, making the biphoton state 
always less entangled then an ideal EPR-pair 
In the case of OAM-entanglement, LG modes with 
different radial indices are in general correlated 
(l8| . It is exactly the Schmidt decomposition pro- 
viding a natural set of modes to study transverse 
entanglement in SPDC, meaning that, apart from 
fundamental interest, it is of interest for high di- 
mensional quantum state engineering. 

Angular Schmidt Modes. Bipho- 
ton state generated in type-I SPDC pro- 
cess has the following form [lj|]: |^) = 

\vac) + J dkidk2^{ki,k 2 )\l) k ^\l) k2 , with k^ 2 
being the wavevectors of scattered photons. In 
the thick-crystal approximation and neglecting 
the walk-off effect caused by birefringence of 
the nonlinear crystal, the biphoton amplitude 
yffeii kg) is described by the following expression 
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where £ p (k) stands for angular spectrum of the 
pump, L is length of the crystal, k p - wavevector 
of the pump, and subscript _L denotes transverse 
vector component. We may rewrite the amplitude 
in a form of Schmidt decomposition: 



*(h±, fc 2 ~l) = v / ^0i(*Tx)V'i(*2x). (2) 

i=0 



Unfortunately, there is no analytical expression 
for Schmidt modes known for an exact wavefunc- 
tion ((T|). We will assume the pump to be gaus- 
sian and approximate the biphoton amplitude by 
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a double-gaussian form as proposed in [23, 24 1 
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which allows one to find Schmidt decomposition 
in an analytical form 0|. The wavefunction © 
depends on two experimentally adjustable param- 
eters: pump divergence a and phase-matching an- 
gular width b, and the corresponding Schmidt 
modes may be shown to be Hermite-Gaussian 
or Laguerre-Gaussian modes with appropriately 
chosen widths. Namely, in cartesian coordinates 
{fa(2)x> fa(2)y} the decomposition reads: 



^ y/ A„ X m Ipn (fax ) 1pm (fay ) X ^ n (fax)lpm(fay) , 
mn 

(4) 



with l/j„(fa x ,2x) = (^) U^fax,2x 

and n (x) = (2™n!0F)- 1 / 2 e- ;r2 / 2 iJ n (a;) being 
Hermite-Gaussian functions. Eigenvalues decrease 
exponentially with growing n: 
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We performed the numerical calculation of eigen- 
values and eigenmodes for exact wavefunction (JTJ) 
for the experimental values of pump and phase- 
matching bandwidths. The results were indeed 
close to those predicted by double-gaussian model, 
validating its applicability in our experimental con- 
ditions. 

Since expression Q is a product of two de- 
compositions depending on fa 2x an d fa,2yj the 
Schmidt number, quantifying entanglement in the 
system, has the form K = K x x K y , with: 
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Experimental realization. Our main goal was 
to realize projective measurements in Schmidt ba- 
sis experimentally. The simplest case is a zeroth- 
order HGoo mode, which can be filtered by appro- 
priate coupling to a fundamental gaussian mode of 
the single-mode fiber. So a single-mode fiber fol- 
lowed by a photon-counting detector would realize 
a projector on the HGoo mode. For higher order 
modes one has to use phase holograms, transform- 
ing the gaussian beam into HG nm mode. This 
transformation may be realized with high efficiency 
using phase only holograms (25- 27 1. When the 
appropriate hologram is chosen, the correspond- 
ing HG mode is transformed into a Gaussian one, 



whereas the others transform to orthogonal ones, 
so only that mode is transmitted through the 
fiber, realizing a desired projection. For Hermite- 
Gaussian modes the holograms look essentially like 
step-wise phase masks introducing a phase shift 
equal to Axg(H n (x)H m (y)) into the beam. 



He-Cd L1 BBO UVM IF L2 BS „ , PM 01 



SMF 




Figure 1: Experimental setup realizing projective mea- 
surements in Schmidt basis (see text for details). Phase 
mask PM is shown as transparent for simplicity, actual 
realization uses double reflection from SLM as shown 
on the inset. 

Schematic of experimental setup is presented in 
FigQ] We used a 325 nm CW He-Cd laser to pump 
a 2mm thick BBO crystal cut for collinear fre- 
quency degenerate Type-I phase-matching. The 
phase-matching bandwidth was estimated to be 
b = 20 mrad. Laser beam was focused into the 
crystal by a 150 mm quartz lens LI, leading to 
divergence a = 5.8 mrad. The particular values 
of parameters were chosen to reduce the Schmidt 
number to a reasonably low value of K ~ 10, so 
that only several modes of lower order were signif- 
icant. Ultraviolet mirror UVM and wide interfer- 
ence filter IF with 40 nm bandwidth were used to 
cut off the pump and additional luminescence. The 
crystal was put in the focal plane of a 140 mm lens 
L2 to collimate the SPDC radiation into a ~ 3 mm 
wide beam which, after being split on a 50/50 non- 
polarizing beam-splitter BS, was coupled with 8X 
microscope objectives 01,2 to single-mode fibers 
SMF placed in objectives' focal planes. The op- 
tical setup was aligned to match the fundamental 
mode of the fiber to a gaussian beam corresponding 
to zeroth-order Schmidt mode for SPDC radiation. 
The fibers were connected to single-photon count- 
ing modules (Prekin- Elmer), followed by a coinci- 
dence circuit. 

We used a reflective Liquid-Crystals-on-Silicon 
phase-only SLM placed in the transmitted arm 
of the BS to display the holograms. The device 
used has 1280x768 square pixels of lOfim size. To 
check the quality of transformation between HG 
modes, we used a 650 nm single-mode fiber cou- 
pled diode laser to produce a gaussian beam, which 
was expanded by a 20X microscope objective to 
match the SPDC mode size. That was a coher- 
ent source similar to HGoo Schmidt mode. The 
holograms were adjusted to minimize the counts 
rate for a transformed beam, ensuring its orthog- 
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onality to the fundamental gaussian mode. We 
obtained minima with visibility greater than 97% 
for modes with n,m < 4. We further checked 
shapes of the modes obtained by scanning the de- 
tecting SMF tip in the focal plane of the cou- 
pling objective. The obtained dependencies of 
counting rate on fiber tip position R(x) showed 
behavior exactly similar to expected convolution 

R ( x ) x I-oc H nm {x/a) exp (-§£J cxp (-^^-j ■ 
Similar behavior observed for SPDC radiation 
when HGqo mode was selected in the reflected 
arm, while one of the orthogonal HG nm modes 
was selected in the transmitted arm. Coincidence 
counts rate R c for holograms with different n and 
m in the transmitted arm are shown in FigJ5Ja). 
Only n = m = case displays high counting rate, 
while all others are suppressed with visibility over 
90%, in good agreement with what is expected for 
Schmidt decomposition. When some higher order 
mode is selected in the reflected arm also, coin- 
cidences are expected to be registered only when 
its counterpart is selected in the transmitted arm, 
which is confirmed by experimental results pre- 
sented in FigI2Jb) . We observe the somewhat lower 
visibility in this case, which may be explained by 
low quality of phase masks used in the reflected 
channel, since phase-step plates made of glass were 
used instead of the SLM in that arm of the setup. 
We believe the anomaly higher counting rates for 
HGqi and HG\\ modes to have their origin in 
the fact that alignment in the vertical plane was 
much more coarse than in the horizontal one, due 
to technical reasons. These results illustrate the 
single-sum property of Schmidt decomposition, i.e. 
demonstrate the pairwise correlations being char- 
acteristic for Schmidt modes. 




Figure 2: Coincidence counting rates for HG m n modes 
selected in the transmitted arm of the setup, while 
HGqo mode (a) and HGio (b) modes are selected in 
the reflected arm. 

When HG nm mode and HGqo modes were se- 
lected in the transmitted and the reflected arms, 
respectively, and the fiber tip was scanned, the 
coincidence rates behave like an expected convo- 
lution of corresponding modes. Moreover, abso- 
lutely similar behavior was observed, irrespectively 
of whether the fiber in the transmitted or in the re- 
flected arm was scanned. This behavior is a direct 
experimental evidence of the fact, that projecting 



one of the photons into a Schmidt mode leads to 
a pure (spatially coherent) state for the second 
one, which is another feature specific for Schmidt 
decomposition. At the same time a very differ- 
ent dependencies were observed for single counts 
of the detector in the transmitted arm - smooth 
wide curves with no local minima, as we would 
indeed expect for spatially multi-mode radiation. 
The value of single counts rate for central posi- 
tions of both fiber tips is directly related to single- 
photon density matrix eigenvalues in Schmidt ba- 
sis: R s ~ (Vw| /0 (s) \ipmn) ~ A m „, providing a di- 
rect way to perform state tomography. Results 
for this kind of measurements are presented in 
FigEJa). These should be compared with eigen- 
values calculated according to the analytical ex- 
pression ([5]) illustrated in Fig^b). To make the 
comparison quantitative, one may use Uhlmann's 
fidelity 

F = Try/ y/pp^y/p = ]T (7) 

m,n 

where Amn are experimental estimates for eigen- 
values obtained by appropriate normalization of 
single counts rates. We achieved a value of F = 
(92 ±3)%, showing a good agreement with double- 
gaussian model. 

Experimental measurement of eigenvalues in 
Schmidt decomposition provides a direct way of 
quantifying entanglement in considered system. 
We have achieved the following experimental es- 
timates for Schmidt numbers: K x = 3.1 ± 0.9 and 
K y = 2.7 ± 0.5 in agreement with theoretical pre- 
diction of K x . y = 2.97 obtained from ((6]). 

The most direct measurement of Schmidt modes 
shape may be performed with ghost-interference 
techniques (28|. For this purpose we substituted 
the single-mode fiber in the reflected channel with 
a multi-mode one with approximately ten times 
larger core diameter. Multi-mode fiber, support- 
ing many spatial modes, served as a "bucket" de- 
tector, collecting the entire angular spectrum of 
SPDC radiation. A 200 /im slit was scanned in 
front of the focusing objective in the same arm to 
obtain an image in the coincidence counts. In this 
case, we can directly resolve the angular shape of 
the desired mode with corresponding hologram in- 
serted in the transmitted arm. Fig[4] shows the 
obtained "ghost" interference patterns. They are 
well approximated by Hermite-Gaussian functions, 
as expected. 

Discussion. We have reported here the first, up 
to our knowledge, experimental attempt to ad- 
dress the physical properties of Schmidt decompo- 
sition for an infinite-dimensional system. We have 
applied the classical technique of spatial modes 
transformation to study transverse entanglement 
in SPDC in Schmidt basis. Although a significant 
amount of work was devoted to orbital angular 
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Figure 3: (a) Single counts rates for HG m n modes se- 
lected in the transmitted arm of the setup. Not nor- 
malized, only background noise was subtracted, (b) 
Single-photon density matrix eigenvalues A mn in HG- 
modes basis, which are the weights in the Schmidt de- 
composition, (c) One-dimensional projection on the 
plane n = 0. Red bars are normalized experimental 
results, black bars - results of numerical calculation 
with exact wavefunction |T]), blue curve - eigenvalues 
for HG-modes predicted by double-gaussian model. 



Besides studying spatial entanglement in SPDC 
itself, Schmidt decomposition provides a natural 
way to get pure spatial states of heralded single 
photons, which is of great importance for quantum 
information tasks. It is thus an important step to 
demonstrate the realization of projective measure- 
ments in Schmidt basis. An interesting theme to 
address is the possibility to prepare an initially fac- 
torized spatial state of biphoton pairs by manipu- 
lating the angular spectrum of the pump. Develop- 
ment of Schmidt mode filtering techniques in fre- 
quency domain would also be an interesting chal- 
lenge. We believe that proof-of-principle results of 
this work to be another step on the way to mas- 
tering high dimensional quantum state engineering 
with spatial states of photons. 

We are grateful to M.V. Fedorov for stimulating 
discussions. 

This work has been funded by Ministry of ed- 
ucation and science of the Russian Federation 
(Minobrnauka), the state contract 02.740.11.0223 
and by the RFBR grants 10-02-00204 and 10-02- 
90036_Bel. S.S.Straupe is grateful to the Dynasty 
Foundation for financial support. 




Figure 4: Ghost images of the first three Schmidt 
modes: HGoo (a), HGio (b) and HG20 (c), corre- 
spondingly. Solid red curves are Hermite-Gaussian fits. 



momentum entanglement of Lagguerre-Gaussian 
modes, the question of conditions under which 
they constitute a Schmidt basis is much less stud- 
ied. Results presented here demonstrate, that 
another basis, namely that of Hermite-Gaussian 
modes, may be equally attractive for applications, 
although it is surprisingly much less used in spa- 
tial entanglement experiments. Cartesian coordi- 
nates on the plane of transverse momentum com- 
ponents of photons, corresponding to HG-modes, 
may be especially convenient in the case of high 
transverse entanglement, when the biphoton am- 
plitude is anisotropic, having different shapes in 
direction of optical axis of the crystal and in or- 
thogonal one [5|. 

Presented results clearly show, that (at least for 
mildly focused pump) a set of Hermite-Gaussian 
modes has all the properties of a Schmidt decom- 
position. Namely, perfect one-to-one correlations 
between modes, spatial coherence and expected be- 
havior of eigenvalues. 
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